Prewetting transitions of Ar and Ne on alkali metal surfaces 
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We have studied by means of Density-Functional calculations the wetting properties of Ar and Ne 
adsorbed on a plane whose adsorption properties simulate the Li and Na surfaces. We use reliable 
ab-initio potentials to model the gas-substrate interactions. Evidence for prewetting transitions is 
found for all the systems investigated and their wetting phase diagrams are calculated. 
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PACS numbers: 68.45. Gd;68.35.Rh 
I. INTRODUCTION 

Direct evidence for the first-order nature of the wetting 
transition, as originally predicted by Cahn 0, together 
with the characteristic prewetting jumps away from the 
coexistence line, has been obtained in several different 
systems: quantum liquid films physisorbed on heavy al- 
kali metals complex organic hquids [|, near-critical 
liquid mercury and binary liquid crystal mixtures 

Quite generally, wetting transitions at temperatures 
above the triple point are expected for weakly attractive 
substrates. In particular, it was argued that alkali metals 
provide the weakest adsorption potentials of any surfaces 
for He atoms and therefore that wetting transitions could 
be observed on such substrates This hypothesis has 
subsequently been confirmed and the corresponding tran- 
sition studied in many laboratories ■ More recently, 
similar phenomena have also been predicted and/or seen 
for H2, Ne, and Hg on various surfaces [pd|-p6[. 

Theory dictates that when a first order transi- 

tion from partial to complete wetting occurs at a tem- 
perature above the triple-point (and below the bulk 
critical temperature), then a locus of first-order surface 
phase transitions must extend away from the vapor-fluid 
coexistence curve, on the vapour side. At T < T^, the 
thickness of the adsorbed liquid film increases continu- 
ously with pressure, but the film remains microscopically 
thin up to the coexistence pressure Psat{T), and becomes 
infinitely (macroscopically) thick just above it. The ad- 
sorption isotherms reach Psat{T) with a finite slope. At 
temperatures <T < {T^ being the prewetting 
critical temperature) the thin film grows as the pressure 
is increased until a transition pressure Ptr{T) < Psat{T) 
is reached. At this pressure a thin film is in equilibrium 
with a thicker one, and a jump in coverage occurs as P in- 
creases through Ptr (T) . At still higher pressures this first 
order transition is followed by a continuous growth of the 
thick film which becomes infinitely thick at Psat (T) . The 
discontinuity (i.e. the difference in thickness between the 
two films in equilibrium) becomes smaller and smaller as 
T increases above , until at a temperature T|'™ the two 
films are no longer distinguishable from one another. Fi- 
nally, when T > T^^ , the adsorbed liquid film increases 
continuously with pressure to become infinitely thick at 



Psat{T). 

Noble-gas fluids (other than He) adsorbed on weakly 
attractive substrates such as the surface of alkali metals, 
are good candidates for showing prewetting transitions. 
There is however a limited experimental knowledge for 
these systems. The measurements of Hess et al. [Q in- 
dicate that Ne undergoes a wetting transition at a tem- 
perature within 2% of Tc on a Rb surface and that it 
undergoes a drying transition on Cs. Non-wetting of Ar 
on a Cs surface has been proposed on the basis of des- 
orption experiments . No experimental information is 
available at present on the wetting properties of Ar and 
Ne on other alkali surfaces. 

In the absence of experimental observations, numerical 
simulations may play an important role in understand- 
ing the wetting properties of fluids and may be used as a 
useful guide in the choice of systems to be studied exper- 
imentally. We present in this paper one such calculation, 
based on a Density Functional approach, where the ad- 
sorption properties of simple classical fluids (namely Ar 
and Ne) on Li and Na surfaces are studied. Ne on Rb has 
also been studied as a test case, where the wetting tran- 
sition has experimentally been observed jl^]. Recently, 
results of extensive Monte Carlo (MC) simulations on the 
Ne/Li system have been presented which seem to 
show that this system is nonwetting for all T between the 
triple point and the critical point. As shown in the fol- 
lowing, our calculations predict quite a different scenario, 
and are in fact consistent with prewetting transitions be- 
low the critical point. We will discuss the discrepancies 
between our findings and the MC results of Ref. in 
Section 3. 



II. COMPUTATIONAL METHOD 

In recent years Density Functional (DF) methods have 
become increasingly popular because of their ability 
to describe inhomogeneous fluids and phase equilibria. 
Comparison with simulation results shows that, once the 
long-range (van der Waals) attractive forces exerted by a 
surface on the gas atoms are included in the free energy 
density functional, the method provides a qualitatively 
(and most often quantitatively) good description of the 
thermodinamics of gas adsorption on a solid surface and 
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in particular it is able to predict correctly a large variety 
of phase transitions (wetting, prewetting, layering, etc.). 

In addition, also the microscopic structure of the fluid 
in the vicinity of a solid surface, e.g. the oscillatory 
behavior of the density profile due to "packing" effects, 
may described accurately by the more sophisticated Den- 
sity Functionals, which treat in a non-local way the 
short-ranged repulsive part of the intermolecular poten- 
tial |2 



-21]. In our calculations we use the DF proposed 
in Ref. [ 26 1 which has proven to describe successfully the 
structure of a classical liquid in very demanding situ- 
ations, e.g. in the presence of hard- walls or strongly 
attractive surfaces, where rapid variations of the liquid 
density occur in the close vicinity of the adsorbing sur- 
face. Although this functional does not seem to describe 
correctly inhomogeneities of infinite spatial extent such 
as those occurring in a freezing transition, it is fairly 
accurate in describing density inhomogeneities of finite 
extent. 

In the Density Functional approach the free energy of 
the fluid is written in terms of the density p{r) of the 
fluid as: 



F[p]^FHs[p] 
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p{r)p{r')ua{\r ~ r'\)drdr' 



+ / p{r)Vs{r)dr. 



(1) 



Here Fhs is the free-energy functional for an inhomo- 
geneous hard-sphere reference system, the second term 
is the usual mean-field approximation for the attractive 
part of the fluid-fluid intermolecular potential Ua [0, 
while Vs [r) is the external adsorption potential due to the 
surface. For Fhs we use the non-local functional of Ref. 
p6| written in terms of a suitable coarse-grained density 
obtained by averaging the true fluid density over an ap- 
propriate local volume. This scheme predicts good triplet 
correlation functions c^^''(r,r') for a bulk one-component 
fluid. In addition, applied to liquid-solid interfaces it 
gives fairly good results, when compared to "exact" 
Monte-Carlo computer simulations, if the nonuniformi- 
ties are not too large |2^. For instance, both DF and MC 
simulations agree in predicting and locating the prewet- 
ting transitions in the case of Ar/C02 system As 
an even more stringent test, the two-dimensional limit 
of the theory of Ref. has been studied : the be- 
havior of a 3D DF theory in this limit is a good test 
of its performances for describing adsorption phenomena 
at low coverages. The overall agreement is surprisingly 
accurate [^ . 

One major weakness of the free-energy functional de- 
scribed above is its mean-field treatment (i.e. its neglect 
of correlation effects) of the long-range attractive part of 
the adatom-adatom interaction appearing in the second 
term of Eq.(|l|). 

This term gives the correct non-local dependence of 
the free energy with respect to variations of the density 
at two points r,r' away from each other, but it neglects 



the correlation effects which should enhance the attrac- 
tive interactions when |r — r'| is close to the minimum of 
the adatom-adatom potential. If one includes this effect 
within perturbation theory, an expression for the free en- 
ergy similar to Eq. (^ should be used, but with the pair 
correlation function of the adsorbed phase multiplying 
the interaction in the integral of the second term. Rather 
than following this approach, which is cumbersome and 
impractical from a computational point of view, we follow 
Ref. and continue to use the much simpler form (|l|) 
but regard Ua as an effective attractive interaction which 
incorporates the main effect of the pair correlation. This 
way of proceeding is found to improve the prediction of 
the DF theory in many cases |2^,^. We follow previ- 
ous prescriptions p9| , ^ , ^ and use in the second term 
in Eq.(j^) an effective interaction: 



Ua{r) —0 , r < Aeo- 



r > AsCT 



(2) 



which for A = I , corresponds to the bare attractive in- 
teraction. For A < 1 this effective potential has the same 
long range properties as the bare interaction, which is 
an important factor when studying wetting phenomena, 
but its value at the minimum is increased by a factor 
A~^, to simulate qualitatively the pair correlation effects 
described above. At any temperature, we treat the HS 
diameter a and the enhancing factor A as free parameters 
determined by requiring that the experimental values of 
the liquid and vapor densities at coexistence, pi and pv, 
are reproduced for the bulk fluid, as explained in detail 
in the following Section. 

As for the interaction with the substrate, we make 
the usual approximation of treating it as an inert, pla- 
nar surface acting on the the fluid as an external po- 
tential Vs{z) {z is the coordinate normal to the surface 

fae). Accurate ab initio potentials are now available 
, which describe the interaction between noble gases 
and the alkali surfaces. Approximating the true surface 
by an ideal plane should not be a very bad approxima- 
tion for these systems, given the uniformity and almost 
complete lack of corrugation of clean alkali surfaces. The 
equilibrium density profile p{z) of the fluid adsorbed on 
the surface is determined by direct minimization of the 
functional in Eq. (|l|) with respect to density variations. 
In practice, at a given temperature, we fix the coverage 
F — J{p{z) — py)dz and solve iteratively the Euler equa- 
tion p = 5F/5p{z) by using a fictitious dynamics, with 
the value of the chemical potential p fixed by F. The 
results of such minimizations are reported and discussed 
in the following Section. 



III. RESULTS AND DISCUSSION 

Our results are based on the analysis of adsorption 
isotherms, calculated by using the DF described in the 
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previous Section, for the case of Ar and Ne adsorption on 
planar surfaces representing Li and Na surfaces, respec- 
tively. 

The bare Ne-Ne and Ar-Ar interactions are usually 
given in the form of a Lennard-Jones (LJ) 12-6 potential, 
with parameters [ ^4|j20|] tNe-Ne/k = 33.9 K, a Ne-Ne = 
2.78 A and €Ar~Ar/k = 119.8 K, a Ar-Ar = 3.4051. 

We use in Eq.(||) the above values for e, while we de- 
termine the two adjustable parameters entering the func- 
tional (|l|), i.e. the HS diameter a and the corrective fac- 
tor to the potential well depth A, by imposing that the 
experimental bulk phase diagram of the adsorbed fluid 
(either Ar or Ne) is correctly reproduced. In particu- 
lar, it is crucial to verify that the choice of these two 
parameters leads to the correct thermodynamic equilib- 
rium conditions. We have calculated several isotherms 
for a bulk system. If the temperature is lower than the 
critical temperature then a blip in the P — V plane, in- 
cluding region of positive slope, {dP/dV)T > , devel- 
ops. As usual, the densities of the coexisting liquid and 
gas and the equilibrium pressure are found by applying a 
Maxwell (equal-area) construction in the P-V plane. For 
a particular choice of a and A, this construction gives the 
vapor pressure P and the liquid and vapor densities p„ 
and pi at a given temperature. At any T, we determine 
the parameters a and A giving the best fit to the exper- 
imental values of P(T), Pv(T) and pi{T)). A summary 
of our best fit parameters are given in Tables I and II, 
while our results for the coexistence line in the {p*,P*) 
plane are shown in Fig. ([l]) (we use reduced units for the 
quantities considered. These are defined as: P* = Pa^/e; 
p* = pa^, where e and a are the LJ parameters of the 
bare atom-atom interaction) . 

We notice that the value of A is such to increase the 
depth of the effective potential Ua over the depth of the 
corresponding bare potential by a factor ~ 1.78 for Ar 
and ~ 1.8 for Ne. These values are of the same order 
of magnitude as that obtained in Ref. by fitting Tc 
for Ar and are reasonable in order to describe the effect 
of the peak in g{r) at the intermediate densities of inter- 
est. A second comment refers to the fitted values of the 
ct's, which turn out to be about 5% smaller than the LJ 
values quoted above, usually employed in Monte-Carlo 
calculations ||2^,^. As a consequence, the coefficients 
of the long range attraction, Cq = Aed^ are reduced by 
about 30% with respect to the corresponding above men- 
tioned LJ values, and are in a much better agreement 
with the most accurate theoretical values As re- 

marked in Ref. however, the LJ parameters have 

nothing of fundamentals and are only determined to pro- 
vide reasonable agreement between "exact" results (say, 
from Montecarlo simulations) and the experimental re- 
sults in the bulk liquids. In other words they are not 
the values which would apply to an isolated pair of Ar 
or Ne atoms. Similarly, our a's are fitted so that the val- 
ues predicted by our DF calculation for the coexistence 
pressure at a given T and for the corresponding densi- 
ties agree with their experimental counterparts. In view 



of this, our effective interaction entering Eq.(^ should 
not be regarded as a pair potential and we do not expect 
that, when inserted in a Montecarlo simulation, it would 
necessarily lead to the same results as those found from 
our DF treatment. 

Once the bulk properties are optimized in the way 
described above, we switch on the adatom-surface po- 
tential in Eq. (|^) and compute the chemical potential 
p corresponding to the equilibrium density profile for a 
given coverage, as described in Section 2. The coverage 
is expressed in nominal layers I = Pi "^^^ [pi^) ~ Pv]dz. 
The collection of pairs (/i, I) at equilibrium, at a given 
T, represents an adsorption isotherm. As an example 
the isotherm of Ar/Li at T = 128K is shown in Fig.(||), 
where the chemical potential, measured with respect to 
its value at coexistence, is plotted as a function of the 
film thickness . For small coverages the chemical po- 
tential increases approaching the saturation value from 
below. For certain temperatures like the one considered 
in Fig.(|^), however, a sudden change in slope is observed 
as the coverage is further increased. After reaching a 
minimum value the chemical potential rises again, gener- 
ating in the p — I plane a van der Waals-like loop. Since 
films for which p has a negative slope are unstable, this 
structure reveals the existence of first-order transitions 
between films of different thickness. The amplitude of 
the transition (i.e. the amplitude of the discontinuous 
jump in coverage) is determined by making an equal-area 
Maxwell construction. The results of such a construction 
are shown in Fig.(H) with a dashed line: two equilibrium 
thicknesses are thus identified, li and I2, and the jump 
between them represents a prewetting transition. 

A few equilibrium density profiles are shown in Fig.(|^), 
for coverages smaller and larger, respectively, than the 
thin-thick film values li and I2 determined above. 

We summarize our results for the Ar/Li system in 
Fig.(^), where the calculated adsorption isotherms for 
this system are shown for a number of temperatures. It 
appears from Fig.(||) that a wetting transition occurs at 
T ~ 124 K, accompanied by first-order prewetting tran- 
sitions at higher T. The amplitude of the prewetting 
transition decreases with temperature, until it vanishes 
at a critical value T^™ ~ 130 K. At temperatures higher 
than a continuous wetting of the surface takes place. 

In the {p — T) plane the occurrence of these first or- 
der wetting transitions is represented by a line leaving 
smoothly the liquid-gas coexistence curve at and end- 
ing at the prewetting critical temperature T^^. On gen- 
eral grounds one expects A/i = p — pcoex ~ —aiT — 
TwY^'^, where the exponent 3/2 is related to the expo- 
nent of the van der Waals tail —C/z^ of the surface po- 
tential. We show in Fig.(|^) our calculated values for A/i, 
together with a fit with the expected analytical form. 

Fig.(^,(0),(||) show the calculated adsorption 
isotherms for Ar/Na, Ne/Li and Ne/Na, respectively. In 
all cases a sequence of prewetting transitions is found, 
below the bulk critical temperature (Tc {Ar) — 150.9iir 
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and Tc{Ne) = MAK). In Table III we summarize the 
wetting and critical pre-wetting temperatures obtained 
from our calculated adsorption isotherms. In parenthesis 
we report the wetting temperature as predicted by Eq. (|^) 
(see the following). 

The last line in Table III refers to the system Ne/Rb 
experimentally investigated by Hess et al. [Q. We cal- 
culated two isotherms for this system, at T ~ 43 K and 
T = UK. We find non- wetting behavior at T ^ 43K, 
while at T = 44 K complete wetting occurs. We thus 
estimate ^ 44 K for this system, but we are not 
able to resolve any prewetting line of transitions between 
these two temperatures. Both the wetting temperature 
and the temperature interval jT^ — r|"^'| ^ 1 K separat- 
ing the nonwetting from the complete continuous wetting 
regimes are in excellent agreement with experiments [ p^ . 
Although the almost perfect agreement between theoret- 
ical and experimental values of should be regarded as 
fortuitous in view of the many approximations contained 
in our calculations ( planar surface, DF treatment, ef- 
fective potential fit etc.), nonetheless it is rewarding to 
find that our treatment correctly predicts the observed 
wetting transition. 

A simple heuristic model has been proposed by Cheng 
et al. where the energy cost of forming a thick film 
is compared with the benefit due to the gas-surface at- 
tractive interaction, resulting in an implicit equation for 
the wetting temperature 

/•OO 

(pi-py) / Vs{z)dz = ~2-f (3) 

J Zo 

Here pi and p^ are the densities of the adsorbate liq- 
uid and vapor at coexistence, 7 is the surface tension of 
the liquid and zq is the equilibrium distance of the gas- 
surface interaction potential Vs . If one uses in Eq. the 
potential K, used in our calculations, and the experimen- 
tal values for pi{T), Pv{T) and 7(r), values of the wetting 
temperature, systematically lower than those found 
in our calculations are predicted, as shown in the column 
of Table III. One issue raised by the present results is 
thus an apparent failure of this simple model to describe 
wetting behavior in the case of ultraweak potentials, such 
as those investigated here. 

Prewetting transitions in classical fiuids are notori- 
ously elusive. On the basis of model calculations, it has 
been suggested that one possible reason for this elu- 
siveness is that they lie so close to the adsorbate bulk co- 
existence line as to render them difficult to detect. This 
happens because the adatom-substrate potential is com- 
parable in strength with the adatom-adatom potential. 
The model calculations of Ref. |^ show that, for a wide 
range of interaction parameters, the prewetting line lies 
at a chemical potential which differs from that of coexis- 
tence by an amount on the order of IQ^^KTc. Our results 
show no exception to this behavior, although the proxim- 
ity to coexistence, which makes prewetting transitions so 
difficult to observe experimentally, is slightly lower than 



expected. In the case of Ne/Li the values of the chemical 
potential for which the discontinuities in coverage are ob- 
served are only IK (at most) below the saturation value 
(see Fig. 7), while for Ar/Li they are 2K (at most) below 
(see Fig. 4). Thus ^ lO^^i^TTc. The vapor pressure 
at which the prewetting jumps should be observed for 
the system we have investigated, may be approximately 
estimated by assuming ideal gas behavior for the vapor 
phase (see Ref. for a better estimate). One finds 
P = P,^texp{Ap/T) ~ 0.985P,at for Afir^2K, thus in 
a range which may be accessible to experiments. 

We mention at this point recent results based on 
extensive Monte Carlo simulations on the behavior of Ne 
on a Li surface (the gas-surface potential used in these 
simulations is the same as used here ) where nonwetting 
behavior up to the critical point has been found, at vari- 
ance with the results presented here. One may invoke 
various explanations for these different findings. For ex- 
ample one could think that metastability in the region 
close to the critical point may affect seriously MC simu- 
lations and that reliable results may require a number of 
MC moves prohibitively large to reach the regime where 
the wetting transition is found. 

Another possible explanation on why Boyan et al. do 
not see wetting of Ne on Li is that the prewetting regime 
is so close to saturation that their simulation cannot dis- 
cern it. In fact at 40 K we find a prewetting jump of 
slightly more than two layers at a chemical potential only 
0.5K below saturation. At higher temperatures, where 
the jump is higher, it occurs even closer to saturation. 
The situation is even worse for Ne/Na. It seems to be 
more promising for Ar/Li, where one should see a ~ 2 — 3 
layers jump 1.5K below saturation at T = 128A', so 
this is the system were we suggest to make experiments. 

We remark at this point the extreme sensitivity of the 
results presented above to the details of the gas-surface 
interaction potential. It has been shown recently [^ , 
for the case of Ar/C02 system, that minor changes in 
the overall shape of the adsorption potential may alter 
dramatically the properties of the adsorbed film. In or- 
der to evaluate the sensitivity of the present results to 
the shape of the substrate potential, we re-calculated the 
wetting diagram of Ar/Li by using a 9-3 Lennard- Jones 
potential to describe the Ar-surface interaction, instead 
than the ab-initio potentials used above [3^ ], and adjust- 
ing the two parameters entering the 9-3 potential in order 
to have the same well depth D and minimum position as 
the ab-initio potential. We found that the wetting tem- 
perature calculated with this interaction shifts upward 
from 123 to 136 K, while the critical wetting temperature 
changes from 130 to 138. This is not surprising, given the 
dependence of the wetting temperature, as clear from the 
implicit definition Eq. (0), from the global shape of the 
potential, and not just from the well depth and position 
of the minimum. 

Of course the wetting properties will also be very sen- 
sitive to the adatom-adatom potential, which determines 
the surface tension explicitly entering Eq. (||). Thus, 
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on these very weak substrates, wetting or non wetting 
above the triple point, but below Tc, is the result of a 
delicate balance between adatom-substrate and adatom- 
adatom interactions. Another possible explanation of the 
discrepancy between our findings and the MC results of 
Bojan et al. may thus be ascribed to the different 
long range behavior of our effective interaction Ua (see 
Eq.(2) and the ensuing discussion) and of the bare LJ 
potential used in the Montecarlo simulations of Ref. . 

Experiments are currently in progress p9| , which will 
hopefully be able to verify the predictions contained in 
this paper. 

ACKNOWLEDGMENTS 

We thank G.Mistura, L.Bruschi, M.W.Cole and 
S.Curtarolo for useful discussions and critical comments. 



[1] J.W.Cahn, J.Chem.Phys. 66, 3667 (1977). 

[2] P. Taborek and J.E. Rutledge, Phys. Rev. Lett. 68, 2184 
(1992); G.B.Hess, M.J.Sabatini and M.H.W.Chan, ibid. 
78, 1793 (1997); D. Ross, P. Taborek and J.E. Rutledge, 
to be published in Phys. Rev. B; G. Mistura, H.C. Lee 
and M.H.W. Chan, J. Low Temp. Phys. 96, 1 (1994). 

[3] H.Kellay, D.Bonn and J.Meunier, Phys. Rev. Lett. 71, 
2607 (1993). 

[4] V.F.Kozhevnikov, D.I.Arnold, S.P.Naurzakov and M.E. 

Fisher, Phys. Rev. Lett. 78, 1735 (1997). 
[5] R.Lucht and C.Bahr, Phys. Rev. Lett. 78, 3487 (1997). 
[6] D.Beysens and D.Esteve, Phys. Rev. Lett. 54, 2123 

(1985). 

[7] E. Cheng, M. W. Cole and W. F. Saam and J. Treiner, 
Phys. Rev. Lett. 67, 1007 (1991); Phys. Rev. B 46, 13967 
(1992); Erratum B 47, 14661 (1993). 

[8] P.-J.Nacher and J.Dupont-Roc, Phys.Rev.Lett. 67, 2966 

(1991) ; P.J.Nacher, B.Demolder, and J.Dupont-Roc, 
Physica B 194, 975 (1994). 

[9] J. Rutledge and P. Taborek, Phys. Rev. Lett. 69, 937 

(1992) ; J.E.Rutledge and P. Taborek, J. Low Temp. Phys. 
95, 405 (1994); Ross, D., Phillips, J. A., and P. Ta- 
borek, J. Low Temp. Phys. 106, 81 (1997); D. Reineh, 
H. Can, S. Herminghaus and P. Leiderer, Czech. J. Phys. 
46, Suppl. SI, 431 (1996). 

[10] K.S.Ketola, S.Wang, and R.B.Hallock, Phys.Rev.Lett. 

68, 201 (1992); R.B.Hallock, J. Low Temp. Phys. 101, 

31 (1995); A.F.G.Wyatt, J.Klier, and P.Stefanyi, Phys. 

Rev. Lett. 74, 1151 (1995). 
[11] E. Cheng, G. Mistura, H. C. Lee, M. H. W. Chan, M. W. 

Cole, C. Carraro, W. F. Saam and F. Toigo, Phys. Rev. 

Lett. 70, 1854 (1993). 



[12] D. Ross, J. E. Rutledge, and P. Taborek, Fluid Phase 
Equihbria 150-151, 599 (1998); D. Ross, J. E. Rutledge, 
and P. Taborek, Science 278, 664 (1997); J. A. Phillips, 
P. Taborek, and J. E. Rutledge, J. Low Temp. Phys. 113, 
829 (1998); J. A. Phillips, D. Ross, P. Taborek, and J. 
E. Rutledge, Phys. Rev. B 58, 3361 (1998); D. Ross, P. 
Taborek, and J. E. Rutledge, Phys. Rev. B 58, R4274 
(1998). 

[13] G. Mistura, H.C.Lee and M.H.W.Chan, J. Low Temp. 

Phys. 96, 221 (1994); D. Ross, P. Taborek, and J. E. 

Rutledge, Phys. Rev. B 58, R4274 (1998). 
[14] G.B.Hess, M.J.Sabatini and M.H.W. Chan, Phys. Rev. 

Lett. 78, 1739 (1997). 
[15] F. Hensel and M. Yao, Eur. J. Sol. Stat. Inorg. Chem. 

34, 861 (1997). 

[16] V. F. Kozhevnikov, D. L Arnold, S. P. Naurzakov, and M. 

E. Fisher, Fluid Phase Equihbria 150-151, 625 (1998). 
[17] H. Nakanishi and M. E. Fisher, Phys. Rev. Lett. 49, 1565 

(1982). 

[18] For reviews, see S. Dietrich, in Phase Transitions and 
Critical Phenomena, edited by C. Domb and J. L. 
Lebowitz (Academic, London, 1988), Vol. 12, p. 1; Liq- 
uids at Interfaces, edited by J. Charvolhn et al. (North- 
Holland, Amsterdam, 1990), Courses 9, 10, and 11. 
(1980). 

[19] W.Friess, T.Brunner and D.Menzel, Surf. Sci. 182 (1994). 
[20] M. J. Bojan, G. Stan, S. Curtarolo, W. A. Steele, M. W. 

Cole, Phys. Rev. E 59, 864 (1999). 
[21] P.Tarazona, Mol. Phys. 52, 81 (1984). 
[22] T.F.Meister and D.M.KroU, Phys. Rev. A 31, 4055 

(1985). 

[23] S.Sokolowski and J.Fisher, Mol. Phys. 68, 647 (1989). 
[24] T.K.Vanderlick, L.E.Scriven and H.T.Davis, J. Chem. 

Phys. 90, 2422 (1989). 
[25] D.M.KroU and B.B. Laird, Phys. Rev. A 42, 4806 (1990). 
[26] E.Kierlik and M.L.Rosinberg, Phys. Rev. A 42, 3382 

(1990) . 

[27] M.P. Allen and D.J. Tildesley, Computer Simulations of 

Liquids, (Clarendon Press, Oxford, 1990), p. 9. 
[28] E.Kierlik and M.L.Rosinberg, Phys. Rev. A 44, 5025 

(1991) . 

[29] Y. Fan and P.A. Monson, J. Chem. Phys. 99, 6897 
(1993). 

[30] E.Kierlik and M.L.Rosinberg, Phys. Rev. A 44 5025 
(1991). 

[31] E.Velasco and P.Tarazona, J.Chem.Phys. 91, 7916 
(1989). 

[32] E.Bruno, C.Caccamo and P.Tarazona, Phys. Rev. A 35, 
1210 (1987). 

[33] A. Chizmeshya, M. W. Cole, and E. Zaremba, J. Low 
Temp. Phys. 110, 677 (1998). 

[34] Hirschfelder, J. O., C. F. Curtiss, and R. B. Bird, Molecu- 
lar Theory of Gases and Liquids, Wiley, New York, 1964. 

[35] K.T. Tang and J.P. Toennies, Z. Phys. D 1, 91 (1986). 

[36] E.Cheng, M.W.Cole, W.F.Saam and J.Treiner, Phys. 
Rev. B 48, 18214 (1993). 

[37] A.K.Sen and C.Ebner, Phys. Rev. B 33,5076 (1986). 

[38] G.Mistura, F. Ancilotto, L.Bruschi and F. Toigo, Phys. 
Rev. Lett. 82, 795 (1999). 

[39] G.Mistura and L.Bruschi, private communication. 



5 



FIG. 1. Phase diagram of bulk-Ar. 
points. Solid line: experimental curve. 



Squares: calculated 



FIG. 2. Chemical potential (measured with respect to its 
value at bulk liquid- vapor coexistence) as a function of cover- 
age for the Ar/Li system at T=128 K. The dashed line shows 
the value of A/i satisfying the equal-area Maxwell construc- 
tion: the thickness of the "thin" and "thick" film in equilib- 
rium at this value are indicated on the x-axis. 



FIG. 3. Density profiles of Ar films on the Li surface, 
plotted as a function of the distance from the surface plane. 
The thin lines show the density profiles at coverages below 
the "thin" film thickness h (see Fig. 2), the thicker lines show 
the density for film thickness larger than the value h shown in 
Fig.2. The origin of the z coordinates is taken at the surface 
plane position. 



FIG. 4. Adsorption isotherms for the Ar/Li system. Afi 
is measured from the saturation value. Squares: T — 124 K; 
Full dots: T = 126 K; Triangles: T = 128, K; Crosses: 
T = 129 K; Stars: T = 130 K; Open dots: T = 132 K. 



FIG. 5. Calculated wetting phase diagram for the Ar/Li 
system. The squares are the calculated points, the dashed 
line is a fit, as described in the text. 



137 0.066 (0.070) 0.089 (0.087) 0.586 (0.586) 0.9438 0.5886 

138 0.069 (0.073) 0.094 (0.092) 0.578 (0.578) 0.9428 0.5897 

139 0.073 (0.077) 0.100 (0.098) 0.569 (0.569) 0.9416 0.5912 

140 0.076 (0.081) 0.106 (0.103) 0.561 (0.560) 0.9404 0.5923 



TABLE II. Same as in Table I, for bulk-Ne. 



T{K) 



P* 



Pv 



Pi 



a /(JNe 



X 



39 0.058 (0.061) 0.074 (0.072) 0.596 (0.596) 0.9537 0.5526 

40 0.067 (0.071) 0.088 (0.086) 0.574 (0.574) 0.9497 0.5536 

41 0.077 (0.083) 0.106 (0.103) 0.549 (0.549) 0.9455 0.5549 

42 0.089 (0.096) 0.127 (0.124) 0.520 (0.520) 0.9413 0.5567 

43 0.102 (0.104) 0.156 (0.138) 0.482 (0.482) 0.9470 0.5503 

44 0.116 (0.123) 0.206 (0.201) 0.417 (0.416) 0.9347 0.5695 



TABLE III. Wetting (T») and critical prewetting (Tp„) 



temperatures from DF calculations, 
reports the prediction of Eq. (fl) 



The column labeled T' 



Gas / Surface 




Tpyj 


rpC 
-'- W 


rjiexpt 


Ar/Li 


123 (0.81 Tc) 


130 


107 




Ar/Na 


136 (0.90 Tc) 


140 


117 




Ne/Li 


38 (0.86 Tc) 


41 


33 




Ne/Na 


40 (0.90 Tc) 


44 


36 




Ne/Rb 


44 (0.99 Te) 




38 


43 



FIG. 6. Adsorption isotherms for the Ar/Na system. 
Squares: T = 137 K; FuU dots: T = 138 K; Triangles: 
T = 139, K; Crosses: T = 140 K. 



FIG. 7. Adsorption isotherms for the Ne/Li system. 
Squares: T = 39 K; Full dots: T = 40 if; Triangles: 
T = 41ii:. 



FIG. 8. Adsorption isotherms for the Ne/Na system. 
Squares: T = 41 K; Full dots: T = 42 if; Triangles: 
T = A3K. 



TABLE I. Values of a and A used to obtain the agreement 
between the experimental and calculated values of coexistence 
densities for bulk-Ar at the saturation pressure, as explained 
in the text. Calculated values are reported in parenthesis. 

T (K) pi p| a/aAr^Ar A 

124 0.036 (0.037) 0.044 (0.044) 0.669 (0.669) 0.9562 0.5810 

126 0.040 (0.041) 0.049 (0.048) 0.658 (0.658) 0.9545 0.5815 

128 0.044 (0.046) 0.055 (0.054) 0.646 (0.646) 0.9520 0.5840 

129 0.046 (0.047) 0.058 (0.057) 0.638 (0.638) 0.9526 0.5838 

130 0.048 (0.050) 0.061 (0.060) 0.634 (0.634) 0.9508 0.5833 
132 0.053 (0.055) 0.068 (0.067) 0.621 (0.621) 0.9485 0.5855 
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